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ABSTRACT 

The origin of rapid spectral variability and certain spectral correlations of the prompt 
gamma-ray burst emission remains an intriguing question. The recently proposed the- 
oretical model of the prompt emission is build upon unique spectral properties of jitter 
radiation — the radiation from small-scale magnetic fields generated at a site of strong 
energy release (e.g., a relativistic collisionless shock in baryonic or pair-dominated ejecta, 
or a reconnection site in a magnetically-dominated outflow). Here we present the re- 
sults of implementation of the model. We show that anisotropy of the jitter radiation 
pattern and relativistic shell kinematics altogether produce effects commonly observed 
in time-resolved spectra of the prompt emission, e.g., the softening of the spectrum 
below the peak energy within individual pulses in the prompt light-curve, the so-called 
"tracking" behavior (correlation of the observed flux with other spectral parameters), 
the emergence of hard, synchrotron-violating spectra at the beginning of individual 
spikes. Several observational predictions of the model are discussed. 

Subject headings: gamma rays: bursts — radiation mechanisms: non-thermal — mag- 
netic fields — shock waves 



Introduction 



Gamma-ray burst (GRB) prompt emission exhibits rapid spectral variability the origin of which 
is not fully understood yet. While some of it — the hard-to-soft evolution at late prompt stages — 



can be attributed to the high-latitude emission in the standard synchrotron shock model (Kumar 



&; Panaitescu 


2000 


Genet &: Granot 


2008 


), the long- acknowledged "tracking" behavior ( 


Bhat, et 



al. 1994 Crider, et al. 1997 Frontera, et al. 2000j ) remains unexplained within this model. The 
tracking behavior is quite intriguing: the low-energy spectral index a follows the observed flux 
(Crider, et al. 1997t [Pothapragada 20051, with harder index corresponding to higher flux, as is 
illustrated in figure [l] (we used data from Kaneko et al.|[2006 ). To our knowledge, the synchrotron 
shock model (e.g., with self-absorption and comptonization) does not seem to generally predict such 
a correlation. The situation becomes even more dramatic when one takes into account that the 
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hardest observed a's cannot be reconciled with the standard synchrotron theory, which constrains 
the lower energy photon index to be a < —2/3 ( Preece et al.|[l998 20001, and that the most of the 
GRB spectra exhibit a ~ — 1 — the value by no means special within the theory (Kaneko et al. 
2006)0 



A recent theoretical model ( Medvedev 2006 1 of the spectral variability of prompt GRBs predicts 



that the above properties naturally follow from the jitter radiation theory (Medvedev 20001 of 
electron emission in small-scale magnetic fields generated in situ. Note that one should not use here 
the results by (Fleishman 2006 1, as we discuss at the end of Section [2j Here we do not limit ourselves 
to a particular model of a GRB: either a kinetically-driven outflow, including baryon-dominated 
and lepton-dominated ones (i.e., with electron-ion and electron-positron pair plasmas, respectively), 
or a magnetic- field-dominated outflow (e.g., driven by Poynting flux), or any combination thereof. 

The prediction that relativistic shocks generate strong magnetic fields via the Weibel-like 



(particle streaming) instability (Medvedev &; Loeb 19991 has been confirmed with a large number 



of numerical PIC simulations in both baryonic and pair plasmas ( 


Silva, et al.|2003 Frederiksen, et al. 


2004 


Nishikawa et al.|2005 


Spitkovsky|2008| Chang et al.|2008 Keshet et al.|2008 


1 . The generated 



fields are random on a very small (sub-Larmor) scale, hence the synchrotron theory is inapplicable 



in such environments while the jitter radiation seems to be validated in these simulations (Hededal 



20051, whose immediate relevance to GRB physics has recently been strengthened (Medvedev Sz 



Spitkovsky 2009 1 . The generation of small-scale magnetic fields has also been discovered in PIC 
simulations of magnetic reconnection in pair plasmas ( Swisdak et al.|2008 Zenitani Hesse||2008 1 . 
We can speculate that a strong energy release, as in GRBs, can produce electron-positron plasmas 
in situ during a reconnection event, even if the plasma is initially lepton-poor. The strong small- 
scale fields are produced via the Weibel instability by the streams of accelerated particles in the 
reconnection exhaust funnels, hence jitter radiation is expected to be produced here. 

In this paper, we study spectral variability and correlations within the jitter radiation frame- 
work using a simple numerical model of an instantaneously illuminated, relativistically expanding 



thin spherical shell segment within a jet ( Pothapragada et al. 2007 Reynolds et al. 20081. We 



interpret the flux-a correlation as a combined effect of temporal variation of the viewing angle, the 
anisotropy of the radiation field in the shell and relativistic aberration. We mostly focus on the 
-Pi/(-E'peak) — 0( correlation because other parameters — the hard index, the peak energy itself -Epeak; 
and the bolometric flux — are more susceptible to the overall GRB energetics, details of electron 
acceleration, etc.; so, unlike a, more factors can affect them and cause mutual correlations. 



'^It is worth mentioning that BATSE time-resolved spectral catalog (Kaneko et al. 20061 is the largest available to 
date and will remain such for at least a few years. Besides, GBM instrument onboard of Fermi-GRO is very similar 
to BATSE in many respects, hence it is natural to compare the results with BATSE data. 
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2. Numerical model 

One usually assumes that the each individual "spike" or "pulse" in the prompt light-curve is a 
single "emission episode" produced by an internal dissipation event (a shock or reconnection) that 
occurs within the ejecta when two thin plasma shells of different speeds and/or field orientations 
collide. The emission region may generally be curved, hence we use the model of the instantaneously 
illuminated thin spherically expanding relativistic shell. Because of the shell curvature, the photons 
emitted along the line of sight, d' = 6 = 0, see Figure [2| arrive at the detector first whereas the 
high- latitude emission, ^' ^ 1, comes later. [Hereafter, 'primed' variables are measured in the shell 
CO- moving (i.e., source) frame and 'unprimed' ones are in the observer's (detector) frame]. Hence, 
an instantaneous 'flash' in the source frame translates to a light curve in the observer's frame (we 
neglect the cosmological redshift effect, assuming z = 0). The angle-dependent spectrum (i.e., 
the volume emissivity, j^, = n'^P{v' ,0'), n'^ being the electron density and P being the ensemble- 
averaged spectral power per electron) of jitter radiation in the source frame, shown in Figure |3j 
thus translates into a time-dependent spectrum in the detector frame. Integration over the emitting 
volume and accounting for time dilation, relativistic Doppler effect and relativistic angle aberration 
yields the energy flux Fi, seen by a distant observer. 

Let us consider emission from a thin segment of a spherical shell expanding radially with the 
Lorentz factor F = (1 — /5^)~^/^, /? = v/c and confined within a jet of opening angle 0jet, its axis 
can be misaligned with the line of sight through an angle ^obs) as shown in Figure [2j The observed 
flux from an optically thin relativistic source has been calculated by Woods Sz Loeb| ( |1999 1 ; Granot 



et al. (19991; Ryde Sz Petrosian (20021 and others. Here we briefly outline the method following 



Woods k Lo"ebl ( p99l ). We use spherical coordinates with a GRB explosion being at the origin; 



the line-of-sight (horizontal axis, ^ = 0) is from the origin toward the detector located some large 
distance r away from it. If the detector plane is orthogonal to the line-of-sight and a ray makes 
an angle a with it, the observed energy flux through a unit surface per frequency is 



j Iu{ct, 4>) cos adfl J d(p J Iy{a,<j))ada, (1) 



where cos a comes from the obliqueness of the ray to the detector plane, which is negligible at 
small angles: cos a ~ 1; also sin a ~ a. Here Iu{ct,(f)) is the intensity along the ray in the 
direction (a, (p) and dfl is the differential solid angle of the emitting region. From simple geometry 
a « rs'mO/d = {r/d){l — fj?)^^'^, where fj. = cos9, so that for given r, ada = {—r'^/d'^)iJ,d^. 

The intensity is calculated as follows, 

diu = juds = jl,'D'^{fi)dr/fj,, (2) 

where the differential path length is ds = {r'^dO'^+dr'^)^^'^ = dr/ fi and we used the Lorentz invariance 
of the quantity ivjv^ ^ so that jv/y"^ = i'^i jv'"^ ■ The observed and emitted frequencies, v and v\ are 
related as 
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time (s) 



• • 
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a 



Fig. 1. — The tracking behavior in BATSE burst 6124. Note that when Fi, at the peak energy 
(top, dark red curve) is high, the photon index a (bottom, blue curve) exceeds the synchrotron 
limit a < —2/3 shown with the vertical line in the left panel. 




Fig. 2. — Cartoon showing the radiating shell surface of the mis-aligned jet. Photons emitted from 
thin annuli or parts thereof centered on the line of sight arrive at a detector simultaneously. 
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Fig. 3. — Co-moving normalized jitter radiation emissivity of an electron averaged over 
the electron ensemble, P{v',6'), at angles (from hard to soft, blue through red) 6' = 
0.001, 0.5, 2, 6, 10, 40, 88 degrees. The frequency is normalized to the jitter peak frequency. 
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with D = [r(l — j3ii)\~^ = r(l + being the relativistic Doppler boosting factor. The emissivity 
in the co-moving frame, j^,(r2^, r, te)-, measured in erg/ (s cm^ Hz sr), is a function of frequency i^', 
position r = {r,9,(j)), time te, and the direction toward the detector measured in the co- moving 
frame ft'^ = {6' , (j)'), where the angles transform due to relativistic aberration as 



cos I 



1 



(4) 



Because of the shell curvature, the photon arrival time, t, is related to the emission time, t^, as 
t = te — rfj,/c, so that t = is for a photon emitted from the origin at te = 0. Combining the above 
equations together, we arrive at a general result analogous to equation (8) of ( [Woods &: Loeb|1999p , 



F,{t) = ^ j d(t> j dfi j r^drV^jl/^{fl'a,r,t + rf,/c), 



(5) 



where limits of integration are set by the geometry of the source, (</>, /^) G jet. In general, F is a 
function of position and time. Our numerical model calculates the spectral flux for an arbitrarily 
shaped source and T{6' , cj)'). For a jet, one gets 9 E [0, 6'jet + 6'obs], (t> £ [<^i(^)) 4'2{G)], see Figure |2j 
The 0-limits are found numerically for each by ensuring to be within the edges of jet cone; in 
the absence of edges G [0, 2tt\. For simplicity, we also assume F = const hereafter. Assuming that 
emission occurs a,t r = R and to = te = R/ Pc and the shell is of infinitesimal thickness, AR 0, we 



take the emissivity as j'^, = n'^P{v', e')e{AR'/2- 



■RWt'e-t'o) ^ KPi^',(^'W-R)S{te-to), 



where Q and 6 are the Heaviside step- function and Dirac (5-function. Here we used that r' = Fr, 
t' = t/T and the identities S{ax) = \a\~^S{x), Q{ax) = sign(o)0(x), and also introduced the 
electron surface density Sg = n'eAR' . Substituting j^, into eq. ^ and integrating over dr one 
obtains the time-dependent (^-function oc 6{fi — (3^^ + tc/R). Upon integration over d/i and d(f) we 
obtain 

F,[t)=F^V\t)P{v'(t),^x'(t))[ct>2{t)-ct>l(t)), (6) 

where Fq, is the normalization constant: Fu{te) = 1 at the peak frequency, T){t) = R/(Tf3ct), 
v'{t) = v/'D{t), fi'{t) = {V{t)T~^ — l)//3, and (j)i,4>2 are found numerically, as is described above. 
For a conical jet, 02 — 4'i have been calculated analytically (loka & Nakamura 2001 Yamazaki et 



al. 2002). 



The jitter comoving spectrum P(i^', 6') depends on the full three-dimensional spatial spectrum 
of the magnetic field and the electron distribution in the emitting shell. It is calculated elsewhere 



in fuh detail, see Medvedev (2000 20061; Reynolds et al. (2008 2009). In order to keep our model 



simple, we assume the field spatial spectra in the shell (e.g., shock plane) and normal to it (i.e., 
along the Weibel current filaments) to be identical and described by a broken power-law peaked at 
~ A^^ — the characteristic field correlation scale. The electron energy distributioEj^ is assumed 



^The electron distribution can be anisotropic, so more emission can be produced at small 6' . At present, no good 
quantitative model of the electron angular distribution is available, so n oc cos(S') can be a good alternative to the 
isotropic model. 
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to be a power- law with the low energy cutoff at 7e, ?^(7) oc 7"^ for 7 > 7e. Hence the peak jitter 
frequency is i^jitter ~ le'^/^B- The resultant spectrum is shown in figure |3] for several angles 6' 
ranging from (head-on shell/shock) to 7r/2 (edge-on shell/shock in the comoving frame), at a 
higher angle 6' > tt/2 the spectrum is identical to that at the angle tt — 6'. The P{i'', 0') spectrum 
above the peak is set either by the electrons, yielding the slope —{p— 1) /2, or by the field spectrum, 
whichever is harder. If it is set by the field (as it is in this study), the variation of the high-energy 
index of the P{v' ^6') spectrum in the range ~2.5-3.5 is real and is associated with the varying 
viewing angle only, neither with the changes of the magnetic field spectrum, nor with the electron 
energy distribution. The slope of the field spectrum at k < ks plays a minor role for the low-energy 
slope of the radiation spectrum. Note that the variation of the lower-energy index between 1 and 
(or and -1 for the photon spectrum) is a benchmark feature of the jitter radiation regime. 

Regarding jitter radiation, the reader we shall be warned on the following. Recently an invalid 



statement has been made (Fleishman 20061 that the spectrum oc below the peak, "valid 
in the presence of ordered small-scale magnetic field fluctuations, does not occur in the general 
case of small-scale random magnetic field fluctuations." It has been demonstrated ( Medvedev|2005 



2006 ) that the statement is flawed since jitter radiation from random magnetic fluctuations with 
a fairly general distribution function (not just "ordered small-scale" fields) does allow for oc 
spectra. Moreover, the entire range from oc to oc is allowed for monoenergetic electrons, and 
even softer spectra can form, depending on the electron energy distribution. Another confusing 



issue is related to the absorption-like oc u spectrum (Fleishman 20061. Such a spectrum is due 
to plasma dispersion and shall be seen in/near the radio band, not in gamma-rays, as the reader 
might incorrectly infer ( Medvedev|2005 1 . We stress once again here that the spectral shape of jitter 



radiation depends, in general, on the viewing angle: it can vary between oc i^^ to oc u^, and this is 
the effect we use in the present study. 

The obtained light-curve Fu{t) represents the full spectrum evolution in time for each individual 
pulse of a prompt GRB and, possibly, an X-ray flare. In order to analyze spectral correlations and 
to compare our model to observational results, we fit the spectrum at each time step, t, with the 
Band function ( Band et al.|1993 1 to obtain the peak energy ii^peak (or, in general, the break energy), 
the low- and high-energy photon indices a and /3spec and the overall normalization A. The fit is 
done in the linear scale of E (not log{E)). This makes the fit more sensitive to the spectral shape 
around Ep^ak and less sensitive to the shape at energies lesser and greater than -Epeak- This is a 
better method, since it is known that in the region E < Ep^ak, the spectral shape is not a simple 
power law. Many GRB spectra were analyzed similarly, hence the comparison with observations is 
facilitated. Since noise is absent in our spectra, in contrast to observations, the values should 
not be compared. In order to mimic the spectral limitations of BATSE, we limit our six-decade 
spectrum to the two-decade spectral window including most of E'peak- values, except for very soft 
spectra at late times. In our correlation analysis, we use the low-energy photon index a and the 
energy flux at the peak energy, F^{Epea.k)i because the latter is a well-defined physical quantity 
characterizing an emission mechanism, in contrast to the flux within a fixed energy band often 
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Fig. 4. — (left panel) — The single pulse light-curve - the time evolution of a (bottom curve, 
left axis) and normalized flux, -Fi/(-£'peak) (top curve, right axis); (right panel) — the correlation 
of Fi,(£'peak) and a. Three cases are plotted: ^obs = (green dots), ^jet/2 (red dots), and Ojet 
(solid line). The feature in F^{Epcak) at i ~ i?/(2cr^) ~ 10 s is due to jitter spectrum feature at 
oblique angles ~ 1; it is a prediction of the model. The overall tracking behavior is evident. The 
synchrotron-violating spectra are seen at early times t < 0.03 s and at very late times t > 1000 s, 
near a jet break. In the right panel, the upper part of the path may be interpreted as the early 
emission phase of an individual emission episode and the vertical segment is the high-latitude 
emission phase. Steepening of the spectrum at very late times is another observational prediction 
of this model. 



Fig. 5. — (left panel) — The binned synthetic energy flux hght-curve with 100 emission episodes 
(red curve, left axis) and the low-energy index a versus time (blue curve, right axis); (right panel) — 
Scatter plot of a vs. normalized flux -Fjy(-E'peak)- The tracking evolution and the positive correlation 
of F,^{Epeak) and a are evident. Most of the points are clustered at a ~ — 1 and there is a smaller 
population of synchrotron- violating points with a > —2/3 — both results are in excellent agreement 
with BATSE spectral data (Kaneko et al. 20061. The inset shows two representative spectral flts 
(red curves) to the model (blue points) at times corresponding to low and high a's; these times are 
marked by vertical dot-dashed lines in the left panel. The flts are accurate everywhere, except the 
high energy edge. 
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used in GRB data analysis. A detailed study of spectral correlations between other parameters 
goes beyond the scope of this paper and will be presented elsewhere (Pothapragada, et al. 2009, in 
preparation). 



3. Results and discussion 

In our model, ^jet and tc = R/{2cT'^) are free parameters; they were chosen to be tc ~ 10 s and 
^jet = 5°; here we also reset t to be the 'time since trigger', i.e., the very first photons arrive at t = 0. 
Figure [4]^left) shows the time evolution of a and -Fi/(-E'peak) within a single pulse (emission episode) 
for three misalignment angles ^obs = 0, ^jct/2, ^jeti the third case appears between the other two 
because we plot the normalized (not absolute) flux. The differences are minimal: the a-curves 
overlap and the F^{Epeak)-cuvves start to deviate around the jet break at ~ lOOtc- The overall 
tracking behavior is evident: the higher flux in the beginning of the pulse corresponds to the larger 
a, i.e., the steeper spectrum below Speak- At early times and at very late times (near the jet break) 
a exceeds the "synchrotron line of death" a = —2/3. This is a prediction of our model. Figure 
|4]^right) represents correlation of the same spectral parameters. The pulse starts off in the upper 
right corner and follows the upper part of the curve to the left — this is the early emission phase. 
At later times t ^ tc, the path turns downward — this is the high-latitude emission phase. At 
very late times, one sees steepening of the spectrum. Because the magnetic field lies predominantly 
in the shell plane (in the internal shock scenario), we suggest that the polarization of gamma-ray 
emission shall increase at times when one sees the shell nearly edge-on in the comoving frame, i.e., 
at the time t ~ = ii/(2cr^). A stronger polarization is expected for a sufficiently misaligned jet, 
~ 1/r < ^jot, or a distorted internal shock front /reconnection layer. The presence of a feature in 
Fi/(£'peak) at t ~ tc is another model prediction, though the exact shape can depend on the actual 
magnetic field distribution. 

Finally, we create a synthetic GRB by overlaying 100 individual pulses (now with tc ~ 1 s) at 
random positions and with random amplitudes drawn from a power-law distribution derived from 
the pulse statistics ( Beloborodov et al. [2000 1. In order to mimic the time-resolved spectral analysis 



done by the BATSE team, we then bin it in time so that fiuence is the same in all bins (the bin 
sizes, of course, vary); the binned energy fiux light-curve is shown in figure [5]^ left) with the red line. 
The bin-averaged spectra are then fit as usual and the dependence of the soft photon index a on 
time is overplotted in the same graph. The inset in the right panel shows that the synthetic spectra 
are well fit with the Band function. Clearly, a tracks fiux in the same way it does in figure [T| The 
scatter plot, figure [5|right), does show a positive F^{Epeak)-0( correlation with some scatter. One 
can also see that there is a small population of spectra with steep (even steeper than synchrotron) 
a's and that many spectra have their soft indices clustered around a = —1, in great agreement 



with observational data (Preece et al. 2000 Kaneko et al. 2006). We stress that this jitter model 



is the only one, to our knowledge, which can explain these facts altogether. 



4. Conclusions 



We developed a model of prompt GRB spectral variability built upon physical understanding 
of the internal structure of a dissipation region. The model presented here is rather crude, because 
it uses very simplified models of the magnetic field spatial spectrum, the electron distribution and 
the simplified relativistic shell kinematics. Yet, the results obtained mimic the observed spectral 
correlations of the BATSE catalog remarkably well without any ad hoc assumptions, thus suggesting 
robustness of the model. We extend our analysis to other spectral parameters and use more 
elaborate models in forthcoming publications. The following predictions can be made. First, one 
shall expect steepening of the spectral index below the peak energy at very late times. Such 
observations are difficult because an afterglow sets in before this and because the peak shifts to 
low frequencies, where (self-) absorption may play a role. Second, a feature in the light curve 
Fy{Epe3ii,t) can be expected around t ~ R/{2cr^). Third, if the jet is misaligned with ^obs ~ l/T 
and/or the shock front /reconnection layer is distorted, one can expect enhancement of polarization 
of gamma-ray emission at times t ~ R/ilcT"^). In magnetic outflows, the ambient field can also 
contribute to polarization and (synchrotron) emission. 
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